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I.  INTRODUCTION 

This  report  discusses  high  frequency  (15-60  kHz)  models  of  surface  forward 
scattering  and  reflection  that  are  appropriate  for  use  in  acoustic  simulations  based 
on  ray  tracing  techniques. 

Typically,  some  of  the  following  properties  of  a  pulse  after  a  single  surface 
interaction  are  needed  for  modeling  purposes:  (1)  the  amplitude  as  a  function  of 
time  in  or  near  the  specular  direction  (for  example,  this  gives  pulse  elongation);  (2) 
the  corresponding  amplitude  statistics;  (3)  the  partitioning  into  scattered 
(incoherent)  and  reflected  (coherent)  components  in  the  specular  direction;  and  (4) 
the  changes  in  the  frequency  spectrum  due  to  the  moving  surface.  At  a  simpler 
level,  (5)  an  effective  reflection  loss  may  be  needed,  e.g.,  for  reverberation  or  one¬ 
way  propagation  simulations.  At  a  more  detailed  level,  (6)  two-point  statistics  at  the 
receiver  may  be  desired.  In  what  follows  an  attempt  will  be  made  to  describe 
models  for  (l)-(5),  but  item  (6)  will  not  be  addressed.  With  regard  to  scattering,  it 
will  be  assumed  that  the  primary  interest  for  simulations  lies  in  representing  the 
scattered  signal  by  an  effective  reflection  loss  due  to  the  difficulty  of  including 
forward  scattering  in  simulations.  Thus,  a  detailed  description  of  the  angular 
distribution  of  the  scattered  intensity  near  the  specular  direction  will  not  be  given. 

The  present  lack  of  comprehensive  experimental  studies  of  surface  forward 
scattering  at  high  frequencies  forces  considerable  reliance  on  heuristic  arguments, 
theoretical  studies,  the  use  of  experimental  results  at  lower  frequencies,  and  the 
use  of  laboratory  experiments  often  at  even  higher  frequencies.  Important 
absorption  effects  of  near  surface  bubbles  due  to  breaking  waves  must  also  be  taken 
into  account. 
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II.  SCATTERING  AND  REFLECTION 

The  terminology  here  will  be  to  restrict  the  term  ’’reflection"  to  the  description 
of  the  coherent*  component  of  the  pressure  field  that  exists  after  a  surface 
interaction.  "Scattering"  is  associated  with  the  incoherent  component.  The 
distinction  between  reflection  and  scattering  is  important  because  the  range 
dependence  of  the  two  components  can  differ  (see  Clay3  for  a  related  discussion). 
The  amplitude  statistics  and  frequency  redistribution  are  also  affected  by  the 
presence  of  a  reflected  component.  Reflection  can  be  described  by  a  coherent 
reflection  loss.  RLC0/l: 


RI*o/i  =  log  R  . 


where  R  is  the  intensity  reflection  coefficient  of  the  surface.  Forward  scattering  can 
be  described6-6  by  a  scattering  strength,  S$  (similar  to  that  used  for 
backscattering): 


S5  =  10  log  cr  , 


where  a  is  the  surface  scattering  cross  section  per  unit  solid  angle  per  unit  surface 
area,  and  is  dimensionless. 

Fortunately,  the  complexities  of  using  a  scattering  strength  or  scattering  cross 
section  can  be  avoided  in  some  cases.  Theoretical  treatments  of  forward  scattering 
using  the  Kirchhoff  approximation  in  the  Fresnel  limit  have  shown7  9  that  when  the 
transmitting  and  receiving  beam  widths  are  sufficiently  wide  (and  the  pulse  length 
sufficiently  long)  to  include  the  entire  "active”  scattering  region  then  the  total 
intensity  after  a  single  surface  bounce,  excluding  the  direct  path  contribution  and 
neglecting  absorption  due  to  bubbles,  is  given  by 


I  total  ~  fcohrront  +  lincoherml 


Q 

(ri  +  Tz2 


(1) 


•  Hie  term  coherence  is  used  here  in  the  usual  sense1'4.  The  coherent  component  of  the  outgoing  acous 
tic  field  has  a  definite  phase  relationship  with  the  ingoing  or  incident  wave;  it  is  given  by  taking  an  average 
over  an  ensemble  of  surfaces  or,  equivalently,  over  an  ensemble  of  pings  The  incoherent  part  fluctuates 
from  ping  to  ping. 
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where 

It  is  the  transmitted  intensity  at  unit  distance, 

7-j  and  r2  are  the  specular  slant  ranges  from  the  transmitter  and  receiver  to 
the  surface, 

A0  is  a  unit  area. 

This  result  can  be  understood  as  a  consequence  of  energy  conservation  when  the 
beam  widths  are  large  compared  to  the  angular  broadening  due  to  scattering. 
Equation  (l)  is  valid  for  large  or  small  surface  roughness  and  indicates  that  for 
many  conditions  the  surface  interaction  can  be  described  by  a  0  dB  loss  for  total 
intensity  when  attenuation  due  to  near  surface  bubbles  can  be  neglected.  A 
discussion  of  the  conditions  when  Eq.  (l)  is  expected  to  apply  will  be  given 
subsequently. 

Note  that  Eq.  (1)  would  be  the  result  for  a  cw  transmitting  source.  For  a  finite 
pulse  width,  the  time  dependence  of  the  received  intensity  requires  specification. 
Both  this  time  dependence  and  the  frequency  spreading  effects  are  conveniently 
described  by  a  scattering  function.10,11  We  will  avoid  these  more  formal 
developments  here  and  simply  write  the  received  intensity  as 

^received )  =  hotali  Oo*  .  (2) 

where 

Jtotai(t)  is  the  total  time  dependent  intensity  (coherent  plus  incoherent) 
following  the  surface  interaction,  neglecting  absorption, 
ab  is  an  absorption  factor  due  to  bubbles. 

Then,  separately,  the  frequency  spreading  due  to  the  surface  interaction  will  be 
specified.  In  Eq.  (2)  we  are  making  the  assumption  that  the  forward  scattered  signal 
is  dominated  by  the  part  that  penetrates  the  bubble  layer,  is  attenuated  in  the 
process,  and  scatters  from  the  surface. 

The  separation  of  the  total  intensity  into  coherent  and  incoherent  intensities  is 
most  meaningful  when  the  size  of  the  Fresnel  zones  on  the  surface  is  large 
compared  to  the  spatial  correlation  length  of  the  surface  roughness.  This  condition 
should  be  satisfied  for  wind  driven  waves.12  In  cases  of  swell  with  no  wind,  however, 
this  condition  may  not  be  satisfied.  Surface  scattering  in  this  latter  case  has  been 
described  in  terms  of  focusing  or  defocusing  from  curved  surfaces,13,14  but  will  not 
be  considered  here. 


APL-UW  7-83 


3 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


ffl.  SURFACE  ROUGHNESS 

The  surface  roughness  parameter*  is  given  by,s 


X  = 


2nh  sin  +  sin  6 z 
A 


where 

h  =  rms  surface  waveheight, 
A  =  acoustic  wavelength, 

=  incident  grazing  angle, 
62  =  outgoing  grazing  angle. 


(3) 


The  roughness  parameter  is  useful  for  classification  of  roughness  regimes.  For  x  < 

0.5,  the  surface  interaction  is  principally  a  reflection;  for  \  ;>  2.0,  scattering 

dominates.  At  intermediate  values  of  x<  the  two  components  can  be  comparable.  By 

taking  x  into  account,  experimental  data  obtained  outside  the  frequency  range  of 

primary  interest  (15-60  kHz)  can  be  used  to  gain  insight  when  corresponding  data 

over  15-60  kHz  are  lacking.  The  rms  waveheight,  h,  which  results  in  a  roughness 

parameter  of  1,  i.e.,  the  transitional  h,  is  shown  in  Figure  1  as  a  function  of 

frequency  for  grazing  angles  of  5°  and  10°  (assuming  8X  -  8Z  and  c  =  1500  m/s).  If  h 

is  large  compared  to  about  0.05  m  (2  in.),  scattering  can  be  considered  dominant. 

The  rms  waveheight  h  can  be  related  to  the  significant  waveheight  H ±  (the  average 

3 

peak-to-trough  height  of  the  largest  waves)  by  H j  =  4.0A,  if  a  Pierson- 

3  3 

Moskowitz16-17  surface  wave  spectrum  is  assumed.  Since  the  open  ocean  typically 
has  larger  significant  waveheights  than  0.2  m  (0.7  ft),  coherent  reflection  tends  not 
to  be  the  important  mode  of  surface  interaction  in  the  open  ocean  at  high 
frequencies.  Coherent  reflection  loss  will  be  considered  in  the  next  section, 
however,  since  fetch  limited  conditions  with  low  surface  roughness  can  be  of 
interest. 


•  The  symbols  g  or  \Tg  are  also  commonly  used  Instead  of  x 
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f  (kHz) 


Figure  1.  The  rms  waveheight,  h,  vs  the  frequency,  f,  which  results  in  a  surface 
roughness  parameter  of  y  =  1  for  specular  grazing  angles,  8,  of  5°  and  10°.  In  the 
Kirchhoff  approximation,  y-  1  yields  an  intensity  reflection  coefficient  of  e~x  (see 
Section  IVj. 
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IV.  COHERENT  REFLECTION  LOSS 

If  the  total  intensity  following  a  surface  bounce  is  divided  into  coherent  and 
incoherent  intensities,  then  the  reflection  loss  RL*,,/,  due  to  the  surface  interaction 
and  associated  with  the  coherent  part  alone  may  be  of  interest.  Whenever  the 
coherent  intensity  is  significant,  bubble  losses  can  be  neglected. 

The  coherent  intensity  reflection  coefficient  given  theoretically  in  the  Kirchhoff 
approximation,  assuming  a  Gaussian  distribution  of  surface  heights,  is1,2,4 

R  =  e-*2  .  (4) 


where  x  is  the  roughness  parameter  with  0t  =  0Z  =  0.  Measurements18  20  have  shown 
Eq.  (4)  to  be  a  good  approximation  for  artificially  constructed  surfaces  when  x  <  1- 
Laboratory  experiments21  have  further  verified  the  single  scatter  Kirchhoff 
prediction  for  reflection  losses  up  to  about  x  —  2.  In  this  work  the  measured  surface 
height  distribution,  which  differed  from  a  Gaussian  form,  was  used  for  obtaining  the 
Kirchhoff  prediction.  The  reflection  loss  was  predicted  correctly  up  to  about  20  dB, 
which  gives  additional  support  for  the  theoretical  assumptions  underlying  Eq.  (4). 
For  large  \  the  reflected  component  becomes  very  small  compared  to  the  scattered 
component  and  difficult  to  measure.  Since  measurements22  have  shown  the 
Gaussian  height  distribution  to  be  a  reasonable  approximation  for  the  sea  surface, 
Eq.  (4)  can  be  expected  to  apply  approximately  to  sea  surface  scattering. 

When  only  the  coherent  or  reflected  intensity  is  being  computed  following  a 
surface  bounce,  rather  than  the  total  intensity,  we  recommend  use  of  Eq.  (4),  which 
leads  to  a  reflection  loss  of 


RLcoA  (dB)  =  305/ 2  A2  sin  2 6  ,  (5) 

with  /  in  kHz,  h  in  m,  and  where  c  =  1500  m/s  has  been  assumed.  The  reflected 
(coherent)  pulse  shows  no  pulse  elongation,  no  frequency  broadening  or  shifting,23 
no  ping-to-ping  fluctuation  (by  definition),  and  is  restricted  to  the  specular  condition 
0t  =  0g. 

An  important  caveat  must  be  placed  on  Eq.  (5)  for  low  grazing  angles.  Direct 
experimental  verification  is  not  available  for  the  very  low  grazing  angles  (0o-5°)  often 
encountered  in  simulations.  Also,  the  theoretical  result  Eq.  (4)  is  known  to  be 
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unreliable  at  low  grazing  angles  because  shadowing  and  multiple  scattering  have  not 
been  included,  and  the  KirchhofT  approximation  breaks  down  as  well. 24,25  Thus,  when 
X  >  0.5,  and  8  <  yc,  where  yc  is  the  rms  surface  slope  angle  (typically  about  8°;  see 
next  section),  Eq.  (5)  has  the  greatest  uncertainty.  For  \  <  0-5,  more  accurate 
perturbation  theory26  30  can  be  used.  Kuperman31  has  computed  reflection  losses  in 
this  way  at  750  Hz  and  1500  Hz  using  the  Pierson-Moskowitz  spectrum,  and  finds 
higher  losses  than  given  by  Eq.  (5).  However,  perturbation  theory  does  show  that 
the  reflection  loss  vanishes  as  8  -*  0  in  agreement  with  Eq.  (5). 
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V.  INCOHERENT  SCATTERING 

We  now  turn  to  a  more  detailed  description  of  the  scattered  intensity.  First 
some  qualitative  features  of  high  frequency  forward  scattering  will  be  introduced 
based  on  heuristic  arguments.  This  will  motivate  an  estimate  of  the  beam  width 
necessary  for  the  use  of  Eq.  (l)  and  its  time  dependent  analogue. 

A.  Implications  of  Beam  Width 

The  angular  distribution  of  the  scattered  intensity  can  be  estimated  using  the 
high  frequency  limit,  i.e.,  by  using  a  "facet  reflection"  model.  Evidence  that  such  a 
model  has  approximate  validity  for  high  frequency  forward  scattering  is  given  in 
Appendix  A. 

In  the  high  frequency  limit,  the  scattered  field  is  determined  by  the  surface 
slope  distribution.  As  discussed  in  Appendix  A.  an  effective  rms  slope  angle  (7 c )  of 
about  fl°  is  representative  of  the  ocean  surface  when  the  wind  speed  U  >  5  m/s. 
This  implies  that  an  incident  plane  wave  with  grazing  angle  8i  >  2-/0  will  have  an  in¬ 
plane*  scattered  intensity  distribution  extending  roughly  over  the  scattered  grazing 
angles  Gs  =  6X  ±  2-y0,  or  typically  with  a  spread  of  about  470  -  30°-35D.  However,  for 
8l  <  7o-  the  angular  extent  of  the  scattered  intensity  will  be  reduced.  For  <  70 
the  scattered  grazing  angle  6S  is  bounded  by  0°  and  roughly  8 *  +  270,  which  implies 
that  the  scattered  vertical  "beam"  width  saturates  at  about  270  ~  16°  as  0t  becomes 
very  small  ( 0*  «  70).  For  70  <  <  270,  the  width  is  roughly  0*  +  270.  These 

estimates  are  consistent  with  AFL  data.32,33 

If  the  widths  of  the  transmitting  and  receiving  beams  are  greater  (ideally  much 
greater)  than  the  angular  width  of  the  scattered  intensity,  then  we  expect  Eq.  (l)  to 
apply.  Thus,  in  typical  propagation  and  reverberation  conditions  where  0*  <  10°, 
vertical  beam  widths  of  200-25°  are  marginally  sufficient.  Multiple  surface  bounces 
will  induce  further  broadening  of  the  scattered  angular  distribution,  and  under 
these  conditions  broader  beam  widths  may  be  needed  for  (1)  to  apply. 

For  vertical  beam  widths  narrow  compared  to  the  angular  width  of  the 
scattered  intensity,  the  intensity  loss  should  exceed  the  0  dB  value  given  by  Eq.  (l), 
but  estimates  of  this  effect  will  not  be  considered. 


•  For  in-plane  scattering,  the  incident  and  scattered  rays  are  in  the  same  vertical  plane. 
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Scattering  out-of-plane,  i.e.,  involving  a  change  in  azimuthal  angle,  is  typically 
confined  to  much  smaller  angular  widths.  For  example,  with  equal  source  and 
receiver  depths  and  equal  incident  and  scattered  grazing  angles,  it  can  be  shown 
that  the  incident  grazing  angle,  0 *,  the  change  in  azimuthal  angle,  <p,  and  the 
corresponding  surface  slope  angle,  y,  are  related  by 

sin  (if/2)  =  tan  0,  tan  y  .  (6) 

If  7  is  set  to  8°  {the  typical  rms  value)  and  if  0,  =  10°,  then  Eq.  (6)  yields  ip  =  2.8°.  In 
this  case,  the  azimuthal  width  of  the  scattered  intensity  would  be  estimated  at  about 
6°.  For  small  grazing  angles  (0,  <  10°),  therefore,  horizontal  (azimuthal) 
transmitting  and  receiving  beam  widths  are  usually  sufficiently  large  so  as  not  to 
affect  the  use  of  Eq.  ( 1). 

Experimental  support  for  Eq.  (1)  is  shown  in  the  results  of  Nichols  and  Senko,34 
who  used  essentially  omnidirectional  sources  and  receivers*  and  found  that  the 
received  intensity  was  independent  of  \  over  the  range  of  0.3  to  about  9.  i.e., 
independent  of  \  from  the  reflection  regime  into  the  scattering  regime.  These 
measurements  were  at  low  frequencies  (400-1500  Hz),  where  absorption  due  to 
bubbles  can  be  neglected,  and  waveheights  were  up  to  9  m  (30  ft). 


•  The  pulse  widths  were  also  long  compared  to  the  pulse  elongation  time  estimates  to  be  discussed  later 
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B.  Effective  Reflection  Loss  for  Forward  Scattering  in  Reverberation  Simulations 

In  reverberation  or  total  backscattering,  ray  paths  involving  surface  bounces 
contribute  to  the  reverberation  level,  particularly  in  shallow  water.  These  surface 
bounces  involve  surface  forward  scattering,  possibly  reflection,  and  absorption  due 
to  near  surface  bubbles.  In  general,  surface  forward  scattering  will  cause  pulse 
elongation  in  the  surface  bounce  signal. 

For  reverberation  simulations,  the  pulse  elongation  due  to  forward  scattering  is 
often  unimportant;  only  the  energy  or  f is  of  significance*.  This  will  be  true 
whenever  the  time  scale  over  which  the  reverberation  level  changes  is  long 
compared  to  the  time  interval  over  which  the  pulse  is  elongated  (see  SectionV.C). 
For  such  a  case,  the  modeling  of  surface  forward  scattering  can  be  simplified 
considerably.  To  see  how  this  occurs,  first  note  that  a  0  dB  intensity  loss  would  be 
implied  by  Eq.  (1)  for  cw  transmission,  assuming  the  beam  widths  are  sufficiently 
broad,  as  discussed  in  the  previous  section.  Added  to  this  would  be  any  attenuation 
due  to  bubbles  (Section  VI).  For  a  very  short  pulse  width,  however,  the  peak 
scattered  intensity  may  not  reach  the  cw  value,  i.e.,  the  important  or  active 
scattering  regions  on  the  surface  may  not  all  contribute  simultaneously.  Even  so,  if 
Eq.  (l)  is  to  hold  for  a  long  pulse,  then  it  follows  that  for  a  short  pulse  there  will  be 
no  energy  loss  at  a  receiver  other  than  spreading  loss  as  in  Eq.  (1)**.  This  justifies 
the  use  of  a  0  dB  effective  reflection  loss  and  no  pulse  elongation  for  modeling 
surface  forward  scattering  in  reverberation  simulations,  whenever  the  time  scale 
over  which  reverberation  changes  is  long  compared  to  the  time  interval  over  which 
the  pulse  is  elongated.  This  should  be  an  adequate  model  for  application  to 
boundary  and  volume  reverberation,  for  example.  The  final  result  is  that  the  surface 
bounce  can  often  be  modeled  with  a  reflection  loss  given  entirely  by  the  "surface 
bubble  loss"  (SBL)  given  by  Eq.  (15)  in  Section  VI. 

These  arguments  have  some  experimental  support.  Adlington35  reported  0  dB 
(±  3  dB)  forward  scattering  energy  losses  for  short  pulse  explosive  charge  data.  In 
this  work  0.4  kHz  £  f  £  6.4  kHz  and  10°  ^  0  £  55°.  No  dependence  on  wind  speed  was 

•  It  la  sufficient  to  take  the  time  interval  of  Integration  here  to  begin  at  the  time  of  specular  path  arrival 
and  to  extend  for  a  pulse  length  plus  several  times  the  characteristic  elongation  time  T ,  which  is  defined 
in  Section  V  C 

**  This  can  be  seen  by  viewing  a  long  pulse  as  a  sequence  of  adjacent  short  pulses.  Superposition  of  the 
scattered  intensities  for  the  sequence  of  short  pulses  will  yield  Eq.  (1)  only  if  each  short  pulse  shows  no  en¬ 
ergy  loss  at  a  receiver  due  to  surface  scattering. 
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noted  over  the  range  of  3-10  m/s  (5-20  kn).  Also,  recent  high  frequency 
reverberation  simulations  at  APL38,27  for  low  sea  state  conditions  have  obtained  the 
best  fit  to  measured  reverberation  using  a  0  dB  effective  surface  reflection  loss 
model.  For  all  these  cases,  no  significant  losses  due  to  bubbles  would  be  expected. 

A  potential  difficulty  with  using  the  0  dB  loss  model  of  surface  scattering  in 
simulations  should  also  be  noted.  In  the  previous  section  it  was  estimated  that  when 
8X  «  7q  the  scattered  intensity  would  be  distributed  over  the  range  of  scattered 
angles  8S  from  0°  to  270  -  16°  for  typical  surface  conditions.  In  other  words,  the 
average  scattered  angle  is  not  equal  to  the  angle  of  incidence.  Usually  in 
simulations  all  the  scattered  energy  would  instead  be  assigned  8S  =  6X,  i.e, 
scattering  would  be  treated  as  specular  reflection.  This  difference  will  cause  the 
simulation  rays  to  have  fewer  surface  and  bottom  bounces  per  unit  range  than  would 
actually  occur,  an  effect  of  particular  significance  in  shallow'  water  environments.  In 
addition,  bottom  reflection  losses  at  very  low  grazing  angles  will  be  used  when  higher 
losses  due  to  the  higher  scattered  angles  are  actually  appropriate.  This  coupling 
between  surface  and  bottom  interactions  was  discussed  by  Marsh  and  Schulkin38  in 
the  context  of  semiempirical  expressions  for  shallow  water  transmission  loss. 
Apparently,  procedures  for  using  the  results  of  Ref.  38  in  ray  based  simulations  have 
not  been  developed.  One  plausible  estimate  on  the  simulation  uncertainty  could  be 
ootained  by  comparing  the  standard  8S  =  0X  result  with  an  alternative  model  with  8S 
constrained  to  y0  whenever  9X  <  yo-  As  discussed  in  Appendix  A,  we  would  estimate 
7o  from  the  wind  speed  U  (in  m/s)  at  a  10  m  height  as  follows: 

tan2  70  =  4.6  x  10"3 ln(  2. If/2 )  .  U  >  1  m/s  (7) 


C.  Time  Dependent  Intensity  After  Forward  Scattering 

For  one-way  propagation  with  surface  bounces,  the  time  dependence  or  pulse 
elongation  of  the  scattered  signal  is  of  interest.  Similarly,  in  applications  such  as 
the  simulation  of  target  returns  with  high  time  resolution  pulses,  it  may  be 
necessary  to  account  for  the  pulse  elongation  effect  of  the  surface  interaction.  This 
aspect  of  surface  scattering  is  difficult  to  treat,  since  it  depends  on  the  geometry  of 
the  scattering  event.  The  theoretical  work  of  McDonald  and  coworkers9,11,39-41  will  be 
taken  as  a  guide  in  formulating  an  expression  for  the  time  dependent  intensity.  In 
general  for  a  single  surface  bounce  path,  both  reflected  and  scattered  signals  will  be 
observed  during  the  first  part  of  the  received  signal  with  a  duration  equal  to  the 
transmitted  pulse  length,  while  subsequently  only  the  scattered  signal  will  be 
present.  The  surface  roughness  will  be  assumed  large  (x  £  2),  however,  since  this 
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limit  covers  most  cases  of  practical  interest  at  high  frequencies.  The  reflected 
component  can  therefore  be  neglected  and  only  incoherent  scattering  must  be 
treated.  The  beam  widths  are  assumed  sufficiently  broad  for  Eq.  (l)  to  apply  for 
cw  transmission,  but  the  topic  here  is  the  buildup  and  decay  of  the  scattered 
intensity  for  finite  pulse  widths.  The  intensity  modeled  will  correspond  to  that 
experimentally  obtained  from  an  average  over  many  pings. 

Consider  first  a  single  surface  bounce  with  incident  and  scattered  slant  ranges 
along  the  specular  path  of  t \  and  r2  as  shown  in  Figure  2a.  Using  the  results  of  Refs. 
11  and  39,  we  find  in  the  high  frequency  limit  (large  roughness)  and  for  Gaussian 
surface  statistics  that  the  time  dependent  intensity  scattered  from  the  surface  with 
the  direct  path  excluded  can  be  written 


Is(t)  =  J H(t ,t')It(t')di'  . 
—  00 


(8) 


where  Is(t)  is  the  time  dependent  analogue  of  I incoherent  in  Eq.  (1),  It(t)  is  the 
transmitted  intensity  at  unit  distance,  and  is  an  intensity  impulse  response 

function  of  the  surface  given  by* 


In  Eq.  (9), 


H(t.t’) 


A0 

(n  +  T-a)2 


@(t-f-tg) 

J. 

[nr(<-r-fs)]8 


e  * 


(9) 


>1  o  is  a  unit  area, 

T  is  a  time  related  to  the  time  scale  of  pulse  elongation  and  is  defined  later, 
ts  is  the  specular  path  travel  time, 

tg  =  (r  j  +  r2)/c,  with  c  the  speed  of  sound, 

0  is  the  unit  step  function,  i.e., 

0(< )  =  0  for  t  <0, 

8(< )  =  1  for  t  >0. 

It  is  straightforward  to  show  from  Eq.  (8)  and  Eq.  (9)  that 


•  Eq.  (9)  must  be  modified  to  apply  at  angles  near  vertical  Incidence. 


> 


APL-UW  7-83 


Figure  Z.  Specular  rays  for  (a)  single  surface  bounce  and.  (b)  surface-bottom-surfac 
bounce  paths. 


APL-UW  7-33 


13 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


fUOdt  = 

—  OB 


■^0  _ 

(r,+rj.)2 


fh(t')di 1 


Thus,  zero  energy  loss  at  the  surface  is  insured  for  all  pulse  lengths.  In  particular, 
for  a  rectangular  transmitted  pulse  of  pulse  length  t0,  i.e.. 


h(t)  = 


t  <  0 

0  <  f  <  To  (10) 

t  >  To 


we  have 


fls(t)dt  =  fls(t)dt  = 


ft  To^o _ 

ir  i  +  r2)2 


(H) 


For  the  rectangular  pulse  of  Eq.  (10),  the  explicit  form  for  the  scattered 
intensity  from  Eq.  (8)  and  Eq.  (9)  is 


Is(t)  = 


ItA  o _ 

(T’l  +  T-a)2 


0 

$[v/^/7r 

4>[  Vt/T~  ]  -$[  y/(T~T0)/T 


tS  0 

OStSto  (12) 

T  S  T0 


where  r  -  t  -  ts,  and  where 


e~‘l  dt 


is  the  probability  integral42  (error  function).  T  can  be  viewed  as  a  characteristic 
elongation  time,  since  /,  (t)  given  by  Eq.  (12)  decays  exponentially  with  a 
characteristic  time  T  when  (t  -  To)  /  T»l.  Near  the  time  of  maximum  intensity, 
however,  the  time  scale  of  intensity  rise  and  fall  is  closer  to  T / 2. 

The  time  T  is  given  by 
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2r,r2  tan27C  J 

1  -e-6"70] 

rj  +  r2  c  l 

J 

(13) 


In  Eq.  (13),  the  first  factor  is  the  expression  for  T  that  follows  from  Refs.  11  and  39, 
where  the  grazing  angle  0  is  implicitly  assumed  to  be  large  compared  to  the  rms 
slope  angle  yc.  Results  of  a  simple  facet  reflection  model  for  arbitrary  0/ja  can  be 
incorporated  to  a  good  approximation  with  the  second  factor,  and  thus  we  adopt  Eq. 
(13)  for  T .  As  an  example,  consider  r\  -  r2,  with  0  =  10°,  70  =  8°,  and  c  =  1500  m/s. 
From  Eq.  (13), 


T  =  9.4r  1  =  4.6R  , 


where  R  is  the  horizontal  range,  T  is  in  ms,  and  and  R  are  in  km.  For  R  =  1  km, 
T  -  4.6  ms;  for  R  =  200  m,  T  =  0.9  ms. 

Almost  all  APL  surface  scattering  data  have  been  obtained  with  narrow  beams 
and  thus  do  not  provide  direct  experimental  tests  of  Eq.  (12).  However,  one  recently 
obtained  data  set43  with  good  signal-to-noise  using  an  omnidirectional  receiver  and  a 
broad  beam  transmitter  is  marginally  suitable  for  comparison.  The  experimental 
conditions  were  as  follows: 

Location:  Dabob  Bay  {—  8  km  west  of  Bangor,  WA) 

Wind  speed:  U  =  4.1  m/s  (8  kn) 

Roughness  parameter:  \  ~  1-22 
Pulse  length:  r0  =  1  ms 
Specular  grazing  angle:  0  =  11.7° 

Horizontal  range:  R=  1.26  km 
Transducer  depths:  114  m  and  146  m 
Vertical  beam  width  of  transmitter:  44°. 

The  large  roughness  condition  is  not  fully  met,  and  this  limits  somewhat  the 
suitability  of  the  data  set.  •  Since  \  >s  not  above  2,  coherent  reflection  will  be  taken 
into  account  by  writing 


( total  (T)  = 


O'!  +  rzf 


0  S  t  §  r0 


T  t  T0 


(14) 


APL-UW  7-83 


15 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


with  Is(t)  given  by  Eq.  (12).  From  Eq.  (7),  y0  2  7.3°,  and  Eq.  (13)  yields  T  =  5.53  ms 
The  time  dependent  intensity  from  Eq.  (14)  is  compared  with  the  result  of  a  100-ping 
average  of  the  data  in  Figure  3.  For  this  comparison,  the  maximum  intensity  from 
Eq.  (14)  was  made  to  coincide  with  the  experimental  maximum  by  appropriate  shifts 
in  the  level  and  time  axes.  For  this  case,  the  time  dependence  of  the  intensity  is  in 
reasonable  agreement  with  the  measurement,  but  further  experimental  verification 
is  needed. 

The  single  surface  interaction  treatment  can  be  generalized  to  multiple 
bounces  following  Tuteur  and  McDonald,11  but  there  is  no  direct  experimental 
support.  Consider  the  surface-bottom-surface  bounce  path  shown  in  Figure  2b.  The 
bottom  bounce  will  be  assumed  a  perfect  reflection  to  simplify  the  discussion.  For  a 
rectangular  transmitted  pulse,  the  time  dependent  intensity  is  given  by  Eq.  (12)  with 


T\  +  r2  "*  T\  +  r2  +  r3  +  r4  , 


and 


T  =  7*,  +  Tz  . 


The  time  T  i  is  found  from  Eq.  (13)  with 


r\Tz  ^  ri(r2  +  r3  +  r4) 
r  i  +  r  g  T\  +  (rg  +  r3  +  r4) 

i.e.,  by  effectively  considering  the  second  surface  bounce  to  be  a  perfect  reflection. 
Similarly,  Tg  is  found  from  Eq.  (13)  with 


r\TZ  (r  i  +  rg  +  r3)r4 

T\  +  rg  (r  i  +  r2  +  r3 )  +  r4 


Generalization  to  any  number  of  bounces  is  straightforward. 

In  terms  of  even  more  approximate  forms,  a  0  dB  effective  reflection  loss  with 
no  pulse  elongation  is  adequate  if  total  energy  is  the  primary  interest.  If  only  peak 
intensity  is  relevant,  an  intensity  reflection  coefficient 
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I 

VI.  ABSORPTION  DUE  TO  NEAR  SURFACE  BUBBLES 


The  significance  of  absorption  due  to  bubbles  in  surface  sound-channel 
propagation  has  been  stressed  by  Schulkin,44  and  more  recently  by  Novarini  and 
|  Bruno,48  but  otherwise  this  aspect  of  surface  forward  scattering  has  received  little 

attention.  Figure  4  shows  in-plane  forward  scattering  strengths33  at  a  wind  speed  of 
20  m/s  where  high  absorption  from  bubbles  is  strongly  suspected.  The  model  curve 
is  based  on  the  high  frequency  limit  (see  Appendix  A  for  a  description  of  the  model 
|  and  for  a  definition  of  the  slope  angle  y)\  we  would  expect  the  curve  to  at  least  be  in 

the  vicinity  of  the  data  if  no  absorption  from  bubbles  took  place.  Reasonable 
agreement  with  this  model  is  found  for  low  wind  speed  as  shown  in  Appendix  A.  We 
take  the  large  difference  between  the  model  prediction  and  the  measurement  to  be 
indirect  evidence  of  high  absorption  due  to  bubbles  in  surface  forward  scattering  at 
a  wind  speed  of  20  m/s. 

Unfortunately,  the  modeling  of  this  attenuation  has  large  uncertainties.  The 
data  base*  on  bubble  populations  must  be  considered  insufficient  for  developing 
empirical  forms  for  bubble  size  and  depth  distributions  as  a  function  of  wind  speed. 

|  In  addition,  studies  of  bubble  distributions  by  Thorpe47  show  complicated  spatial 

structure,  which  indicates  that  the  relationship  between  the  average  bubble 
distributions  and  absorption  in  forward  scattering  may  not  be  straightforward. 
Therefore,  the  "surface  bubble  loss"  model  developed  here  has  the  potential  for 
t  large  error.  The  surface  bubble  loss  (SBL)  is  related  to  in  Eq.  (2)  by  SBL  (dB)  = 

'  -10  log  a6. 


The  choices  made  in  arriving  at  a  bubble  loss  model  are  outlined  in  Appendix  B. 
The  loss  estimates  are  based  on  bubble  distributions  measured  by  Johnson  and 
Cooke48  and  Kolovayev.49  The  final  result  is 


SBL(dB)  = 


1.12  x  IQ-7 
sin  8 


U66f 


0.5 


15kHz  g  /  ^  24kHz  (15a) 


I 


SBL(dB)  = 


9.53  x  IQ-10 
sin  8 


U 65  /  20  , 


24kHz  g 48kHz  (15b) 


for  grazing  angle  8,  wind  speed  U  in  m/s  at  a  10  m  height,  and  frequency  /  in  kHz. 

I  _ 


•  Bubble  measurement  techniques  and  results  have  been  summarized  by  Wu  49 
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Figure  4.  Measured  inplane  surface  forward  scattering  strengths  at  high  wind  speed 
compared  with  a  model  based  on  the  high  frequency  limit.  See  Appendix  A  for  a 
description  of  the  model  and  for  a  definition  of  the  slope  angle  y. 


APL-UW  7-B3 


19 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


The  bubble  data  do  not  support  a  sharp  break  in  the  frequency  dependence  at  24 
kHz,  but  there  does  appear  to  be  a  gradual  change  to  a  slower  frequency 
dependence  at  lower  frequencies. 


It  is  interesting  to  compare  Eq.  (15a)  to  AMOS  results44  for  intensity  loss,  w'hich 
Schulkin  suggested  were  due  at  least  in  part  to  bubble  absorption.  Starting  w'ith  Eq. 

(15a),  we  use  the  Pierson-Moskowitz  spectrum17  to  convert  wind  speed  to  rms 
r  !| 

waveheight  h  in  m  U  =  17.7 hz  ;  &  is  estimated50  to  be  2a-3°,  so  we  set  6  =  2.5°;  and 
converting  to  the  units  of  Ref.  44  we  find 


SBL(dB)  =  1.31A2-75(/i/)0-5  . 


h  in  ft 
/  in  kHz 


while  AMOS  data  imply 

Loss  (dB)  =  1.64(A/)0  5  . 


h  in  ft 
/  in  kHz 


The  w’aveheight  dependence  is  very  different,  but  at  rms  waveheights  of  about  1  ft 
(w'hich  corresponds  to  reasonably  typical  wind  speeds  of  about  8  m/s  -  15  kn  ),  the 
losses  are  similar  in  magnitude. 


The  surface  bubble  loss  is  a  very  rapid  function  of  wind  speed  in  the  model  given 
by  Eq.  (15).  Figure  5  shows  the  SBL  just  at  the  onset  of  significant  loss.  An 
appropriate  surface  scattering  data  base,  however,  is  simply  not  available  at  present 
for  verifying  the  validity  of  these  predictions.  For  wind  speeds  above  12  m/s  the  SBL 
predicted  by  Eq.  (15)  for  low  grazing  angles  ( 9  <  10°)  becomes  infinite  for  all 
practical  purposes.  A  low  level  scattered  signal  from  the  bubble  layer  will  remain, 
but  it  is  difficult  to  incorporate  this  effect  correctly  in  an  attenuated  effective 
reflection  loss  model.  A  rough  estimate  based  on  scattering  from  bubbles  (see 
Appendix  B)  indicates  that  the  absorption  loss  would  be  bounded  at  about  30  dB  due 
to  scattering.  However,  such  an  estimate  is  uncertain  because  of  dependence  on  the 
pulse  length  and  beam  widths.  Until  better  information  becomes  available,  we 
recommend  simply  using  a  maximum  attenuation  of  30  dB  whenever  Eq.  (15)  w'ould 
give  an  SBL  >  30  dB. 


Present  modeling  of  surface  bubble  losses  must  be  considered  very  uncertain. 
The  possibility  of  wind  gusts  and  differences  in  the  mean  wind  speed  variation  with 
height  near  the  surface  produce  uncertainties  in  relating  the  wind  speed  at  10  m  to 
bubble  distributions.  Further,  the  likely  occurrence  of  horizontal  variations  in  the 
near  surface  bubble  distributions  will  yield  uncertainties  in  relating  average  bubble 
distributions  to  absorption  in  surface  forward  scattering.  The  surface  bubble  loss  in 
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the  6-8  m/s  wind  speed  range  is  especially  uncertain,  because  this  is  near  the  onset 
of  significant  loss  and  the  predicted  variation  with  wind  speed  is  rapid. 

One  further  point  should  be  noted.  Because  of  their  rapid  falloff  with  depth,  the 
bubble  distributions  measured  by  Johnson  and  Cooke  and  by  Kolovayev  were  almost 
certainly  produced  by  breaking  waves  at  the  surface.  Medwin51,52  has  measured 
bubble  distributions  in  Monterey  Bay,  which  apparently  are  of  biological  origin  and 
thus  are  not  directly  related  to  the  wind  speed.  These  bubbles  are  of  sufficient 
density  to  be  an  important  absorption  mechanism  also,  but  are  not  included  in  the 
model  discussed  here. 
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VII.  FREQUENCY  SPREADING 

The  spectral  properties  of  surface  forward  scattering  have  generally  been 
studied  with  low  to  moderate  surface  roughness  at  lower  acoustic  frequencies  than 
of  direct  interest  here.  From  this  experimental  work53  60  and  other  theoretical 
analyses,61 68  the  following  picture  has  emerged:  At  low  roughness  (x  <  1)  the 
surface  bounce  signal*  is  made  up  of  a  reflected  component,  with  the  same 
frequency  content  as  the  incident  pulse,  and  of  a  scattered  component  forming 
frequency  sidebands  above  and  below  the  incident  center  frequency.  The  sidebands 
are  linearly  related  to  the  surface  wave  spectrum  and  effectively  map  out  this 
spectrum.  The  sideband  maxima  are  displaced  from  the  original  center  frequency 
by  approximately  the  frequency  of  the  peak  of  the  surface  wave  spectrum,  and  the 
magnitude  of  the  sidebands  becomes  comparable  to  the  reflected  component  as  \ 
approaches  1.  (See  in  particular  the  work  of  Roderick  and  Cron.53) 

As  the  surface  roughness  increases  beyond  1,  the  scattered  spectrum  evolves 
into  a  Gaussian  form,  the  width  becoming  large  compared  to  the  frequency  of  the 
surface  wave  spectral  peak,  but  still  small  compared  to  the  corresponding  spectral 
width  for  backscattering.  McDonald  and  Schultheiss83  have  shown  with  theoretical 
studies  that  the  asymptotic  Gaussian  form  is  a  good  approximation  to  the  scattered 
frequency  spectrum  when  x  reaches  2.5. 

Swarts  and  Eggen65  have  computed**  the  width  of  the  Gaussian  frequency 
distribution  as  a  function  of  wind  speed  U  assuming  a  Pierson-Moskowitz  surface 
wave  spectrum.  Converting  their  result  to  the  width  of  the  half  power  points  A/3dB 
leads  to 

A/3dS(Hz)  =  0.1281//  sin  6  ,  (16a) 


•  The  incident  pulse  is  assumed  here  to  be  a  narrowband  signal 

••  The  authors  used  the  KirchhofT  approximation  and  also  invoke  the  Fraunhofer  approximation 
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with  U  in  m/s’  and  /  in  kHz,  or  equivalently, 


A/3dB(Hz)  =  1.76h 2 /sin  6  , 

i 


(16b) 


with  h  in  m.  In  this  model  the  frequency  spectrum  is  broadened  but  unshifted. 


Roderick70  measured  the  forward  scattered  frequency  distribution  for  6  -  10° 
at  20  kHz  and  30  kHz,  and  at  rms  wave  heights  of  0.07  m  and  0.24  m.  The  frequency 
distribution  was  found  to  be  approximately  Gaussian  with  no  evident  shift.  Table  1 
gives  the  3  dB  widths  from  these  data  and  the  corresponding  values  from  Eq.  ( 1 6b); 
the  agreement  is  satisfactory.  For  the  case  with  h  =  0.07  m,  the  measured  surface 
wave  spectrum  showed  significant  deviations  from  the  Pierson-Moskowitz  spectral 
form,  and  some  differences  from  the  predictions  of  Eq.  (16)  would  be  expected  from 
this  alone. 


Table  I.  Comparison  of  Experimental  and  Model  Spectral  Widths 


/  (kHz) 

h  (m) 

Experimental8  3  dB 
Width  (Hz) 

Modelb  3  dB 
Width  (Hz) 

20 

0.07 

1.1 

1.6 

20 

0.24 

2.7 

3.0 

;  30 

0.07 

1.4 

«  , 

30 

0.24 

L _ ™ _ 

4.5 

■  Frorr.  Ref.  70. 
bFromEq.  (16b) 


•  For  direct  use  of  the  Pierson  Moskowitz  surface  wave  spectrurr..  the  wind  speed  should  be  measured  at  a 
height  of  19.5  m.18  However,  use  of  wind  speeds  at  the  standard  height  of  10  m  for  V  in  Eq  (16a)  should 
be  sufficiently  accurate  Wind  speed  normalization  to  a  standard  height  is  discussed  in  Ref  69 
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Eq.  (16)  has  not  been  experimentally  verified  over  a  broad  range  of  conditions 
at  high  frequencies  —  in  particular,  not  at  very  low  grazing  angles  where  Eq.  (16)  is 
known  to  be  unreliable.65  The  reasons  for  uncertainty  are  the  same  as  mentioned  in 
Section  IV:  shadowing,  multiple  scattering,  and  the  breakdown  of  the  KirchhofT 
approximation  at  low  grazing  angles.  Until  such  time  as  better  information  exists 
for  very  low  grazing  angles,  we  recommend  as  an  alternate  model  substituting  the 
effective  rms  surface  slope  angle,  70,  for  ®  in  Eq.  (16)  when  8  <  70-  Some  motivation 
of  this  choice  is  given  in  Section  V.B.  The  slope  angle  70  is  given  by  Eq.  (7). 

Our  recommendation  for  forward  scattering  frequency  broadening  is  to  assume 
an  unshifted  Gaussian  spectral  form  for  all  wind  speeds,  since  conditions  of  low- 
roughness  (x  <  1)  are  relatively  rare  at  high  acoustic  frequencies.  The  3  dB  spectral 
width  is  given  by  Eq.  (16).  but  the  alternate  model  at  low-  grazing  angles  as  noted 
above  is  suggested  as  a  bound  on  the  uncertainty. 

At  high  sea  states  the  effects  of  near  surface  bubbles  must  be  taken  into 
account  as  discussed  earlier.  The  dominant  effect  on  forward  scattering  will  be 
attenuation,  but  scattering  from  bubbles  will  also  occur  and  this  process  may  affect 
the  frequency  spectrum  when  the  attenuation  is  high.  McDaniel  and  Gorman71  have 
recently  considered  the  effect  of  bubbles  on  the  backscattered  frequency  spectrum, 
but  the  forward  scattering  problem  has  yet  to  be  treated.  No  quantitative  estimate 
of  this  effect  is  presently  available. 
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Vm.  STATISTICS 

Several  investigators  have  shown  that  when  the  roughness  is  large  (x  >  2)  the 
amplitude  statistics  over  an  ensemble  of  pings  after  a  forward  bounce  are  described 
to  a  reasonable  approximation  by  a  Rayleigh  distribution. 12-14-72'74  For  low  roughness 
(x  <  2).  the  constant  reflected  component  plus  a  random  scattered  component  has 
been  modeled  with  the  Rician  distribution.12,76-76  Since  the  low  roughness  condition 
seldom  occurs  at  high  frequencies,  and  since  the  Rician  distribution  is  more 
concentrated  about  the  mean*  than  the  Rayleigh  distribution,  we  recommend  use  of 
the  Rayleigh  distribution.  Thus,  the  relative  probability  of  obtaining  a  pressure 
amplitude  (or  envelope  value),  A,  is  given  by77 


P(A)  = 


exP  \-az/2<pz> 


(17a) 


where  p  is  the  instantaneous  pressure  and  <  >  denotes  ensemble  average  at  a  given 
time  in  the  ping  cycle.  The  mean  square  instantaneous  pressure  <pz>  = 

The  equivalent  intensity  distribution  is 


P(I)  =  7j-  exp  [-//</> 


(17b) 


where  /  is  proportional  to  A2.  It  should  be  understood  that  the  models  in  Section  V 
pertain  to  the  mean  </  >. 


•  The  point  is  that  the  Rayleigh  distribution  gives  a  greater  probability  of  large  deviations  frorr,  the  mean, 
and  presumably  it  is  the  more  conservative  choice  to  assume  a  Rayleigh  distribution  for  all  conditions 
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APPENDIX  A 


A  FORWARD  SCATTERING  MODEL 
IN  THE  HIGH  FREQUENCY  LIMIT 


A  simplified  model  for  forward  scattering  strengths  can  be  obtained  in  the  high 
frequency  limit,  i.e.,  by  letting  X  -»  0.  For  a  Gaussian,  isotropic  distribution  of 
surface  slopes,  the  scattering  strength  can  be  written8, 78,79 


=  10  log 


sec4  y 

tan8  y 

4ns 8  6XP 

s8 

. 

(Al) 


where  tan  y  is  the  necessary  surface  slope  at  the  local  scattering  region  that  would 
"reflect"  a  ray  from  the  transmitter  to  the  receiver.  The  mean  square  surface  slope* 
is  given  by 


s2  =  < 


K 

dx 


dy 


>  =  tan8  70  . 


where  £  (x,y)  represents  the  surface  height  and  where  x  and  y  are  orthogonal 
coordinates  in  the  plane  of  the  mean  surface;  the  angle  yo  is  the  rms  slope  angle. 

To  estimate  s8  we  use  the  results  of  optical  glitter  measurements  of  ocean 
surface  slopes  made  by  Cox  and  Mu'.k80,81  and  conveniently  presented  by  Phillips.88 
We  expect  the  high  surface  wave  numbers  to  have  little  effect,  because  of  the  finite 
acoustic  wavelength,  and  therefore  not  contribute  to  the  effective  rms  slope.83  85 
This  is  crudely  accounted  for  here  by  using  the  slope  distribution  measured  with 
surface  oil  slicks,  which  eliminated  surface  wavelengths  shorter  than  ~  30  cm.  The 
same  approximation  has  been  found  useful  in  modeling  electromagnetic  scattering 
from  the  ocean  surface.86  From  Ref.  B2  we  then  obtain  for  wind  speed  U  in  m/s  at  a 


•  Care  is  necessary  here  since  some  authors  define  s2  to  be  the  mean  square  slope  for  one  direction  only 
In  this  case,  sz  in  Eq.  (Al)  and  Eq.  (A2)  is  replaced  by  2sz  =  2  tan 


V 
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10  m  height 


s2  =  4.6  x  10_3ln  (2.l£/2)  ,  (A2) 

for  U  >  1  m/s.  From  Eq.  (A2)  it  can  be  readily  shown  that  the  rms  slope  angle 
saturates  at  about  8°  for  U  £  5  m/s  (see  also  Figure  4.23  in  Ref.  82).  Thus  yo  -  8D  is 
an  estimate  for  an  effective  rms  surface  slope  for  typical  ocean  conditions. 

In  Figure  A-l,  the  scattering  strength  model  given  by  Eq.  (Al)  and  Eq.  (A2)  is 
compared  with  forward  scattering  data  obtained  by  APL.43  Since  for  these  data  the 
scattering  was  in-plane,  i.e.,  the  scattering  region  was  in  the  vertical  plane 
containing  the  transmitter  and  a  highly  directional  receiver,  the  slope  angle  is  given 
6S  ~  6i 

simply  by  y  - - — - ,  where  8S  and  6i  are  the  scattered  and  incident  grazing 

angles,  respectively.  For  specular  scattering,  y  =  0. 

The  approximate  agreement  between  experiment  and  model  in  Figure  A-l  (and 
other  similar  comparisons,  including  APL  vertical  incidence  data43)  supports  the  use 
of  the  high  frequency  limit  for  making  estimates  of  forward  scattering  properties. 


APL-UW  7-83 


27 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


Figure  A-1 .  Measured  in-plane  surface  forward  scattering  strengths  at  low  wind 
speed  compared  with  model  based  on  the  high  frequency  limit.  The  surface  slope 
angle  y  is  zero  for  scattering  in  the  specular  direction. 
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APPENDIX  B 

A  BUBBLE  LOSS  MODEL 


Suppose  an  incident  plane  wave  of  intensity,  Iit  interacts  with  a  surface  with 
incident  and  outgoing  grazing  angles  equal  to  6.  Resonant  bubbles,  assumed  here  to 
be  homogeneous  in  the  horizontal  plane,  will  by  themselves  reduce  the  final 
intensity,  If,  to 

OB 

If  =  /t  exp  -~^fse{z)dz  ,  (Bl) 

where  is  the  total  extinction  cross  section  per  unit  volume  (absorption  plus 
scattering).  Obviously,  near  surface  bubbles  are  not  just  dependent  on  the  depth,  z, 
but  sufficient  information  is  not  available  to  model  the  horizontal  structure.  Note 
that  in  Eq.  (Bl),  the  acoustic  path  is  integrated  up  to  a  depth  d  below  the  surface, 
and  not  to  the  mean  surface  (d  =  0).  For  high  bubble  attenuation,  some  account 
should  be  made  for  scattering  from  wave  troughs,  which  extend  down  from  the  mean 
surface.  This  will  be  done  later  (very  approximately)  by  choosing  a  d  >  0. 

We  define  a  surface  bubble  loss  (SBL)  to  be 

SBL(dB)  =  10  log  =  “10  log  a.b 

(see  Eq.  (2)),  and  Eq.  (Bl)  yields 

<B2> 

The  total  extinction  cross  section  per  unit  volume  is  given  by*2 

S,(z)=  1.7x2Tr2a^ip  /6r  .  (B3) 
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The  factor  1.7  is  given  by  Medwin82  to  account  for  nonresonant  contributions  to 
absorption  and  scattering.  The  radius  of  resonant  bubbles  aR  in  /zm  at  depth  z  in  m 
is  related  to  frequency  /  in  kHz  by 

i 

aR  =  3.25  X  103(1  +  O.lz)2  //  (B4) 


The  quantity  j  is  the  number  of  l  bles  per  unit  volume  per  unit  bubble  radius 

evaluated  at  the  radius  aR,  and  is  discussed  below.  The  damping  constant  6R  is 
given  by 


6r  -  0.0136(1  +  O.lz)  .  (B5) 

Until  better  information  is  available,  we  will  estimate  the  bubble  distribution 
from  the  data  of  Johnson  and  Cooke48  and  Kolovayev49  in  a  form  similar  to  that  used 
by  Novorini  and  Bruno.48  Thus, 


dn 

da 


S0a-> 


(B6) 


with  bubble  radius  a  in  fun.  Below  a  wind  speed  U  of  about  3  m/s,  bubble 

generation  would  not  be  expected  at  all,87  but  for  U  <3  m/s  the  SBL  derived  here  is 
insignificant. 

Johnson  and  Cooke  data  imply  p  =  5*  for 

68  /zm  £  a  £  136  /zm 
24  kHz  <  /  <  48  kHz 

where  here  (and  in  what  follows)  the  frequency  range  is  the  range  of  resonant 
frequencies  that  corresponds  to  the  range  of  bubble  radii.  In  the  Johnson  and  Cooke 


•  This  is  outside  the  range  for  which  Medwin's  factor  of  1  7  in  Eq  (B3)  is  thought  to  be  valid  This 
point  deserves  further  investigation. 
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data  there  is  an  indication  of  smaller  p  for  a  >  136  pm  (/  <  24  kHz).  Kolovayev, 
with  more  data  on  larger  bubbles,  finds  p  -  3.5  for 

1 15  pm  <  a  <  350  pm 
9  kHz  <  /  <  20  kHz 

fv  «  rv 


In  Eq.  (B6)  we  choose 


68  pm  g  a  §  136  pm  (24  kHz  =  f  t*  48  kHz) 
136  pm  i  a  &  280  pm (15  kHz  S  /  S  24  kHz) 


Next,  following  Wu,46  we  take  in  Eq.  (B6)  q  =  4.5,  which  gives  the  dependence  of 
the  bubble  density  on  wind  speed.  This  dependence  is  not  well  determined 
experimentally. 

The  exponential  depth  dependence  of  the  bubble  density  (to  at  least  a  3  m 
depth)  is  consistent  with  the  data  of  Johnson  and  Cooke,48  Kolovayev.49  and  Thorpe.47 
For  the  Johnson  and  Cooke  data,  L  -  1  m  for  U  -  12  m/s.  We  assume  L  is 
proportional  to  the  waveheight.  Using  the  fact  that  for  a  Pierson-Moskowitz  surface 
wave  spectrum  the  waveheight  scales  with  Uz,  we  obtain 

L=(U/  12)z  .  (B7) 


with  L  in  m.  With  these  choices  in  Eq.  (B6)  and  using  the  data  of  Johnson  and  Cooke, 
we  find  S0  =  2.6  x  1010m~3/zm_t  for  p  =  5,  and  Sn  =  1.6  x  107m~3//m_1  for  p  =  3.5.* 

We  choose  d  in  Eq.  (B2)  to  be  given  by  L  to  account  roughly  for  the  presence  of 
wave  troughs.  With  this  choice  d  is  about  1.4  times  the  rms  waveheight  (h),  but  still 
less  than  half  the  significant  waveheight  of  about  4h  (peak-to-trough).  Such  a 
normalization  should  clearly  be  done  by  comparison  with  data,  but  the  required  data 
base  is  lacking. 


•  The  wind  speeds  for  the  Johnson  end  Cooke  data  were  not  measured  at  the  standard  height  of  10  m,  and 
to  be  rigorous,  wind  speed  normalization  to  10  m  should  be  made  (see  Ref.  69).  However,  very  preliminary 
comparisons  with  new  APL  data  (4/83)  appear  to  support  Eq.  (15)  as  is,  taking  V  to  be  the  wind  speed  at 
10  m 
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With  the  rapid  decrease  in  babble  density  with  depth  implied  by  Eq.  (B6)  and 
Eq.  (B7),  the  factor  (1  +  0.  lz)  can  be  set  to  1  in  Eq.  (B4)  and  Eq.  (B5).  Finally, 
combining  Eqs.  (B2)-(B7)  and  doing  the  elementary  integration  over  z  gives  Eq.  (15). 

Scattering  from  bubbles  will  effectively  place  a  limit  on  the  absorption  loss,  and 
this  limit  can  be  estimated  from  the  scattering  properties  of  bubbles.  In  this  limit, 
the  bubble  density  is  so  high  that  the  surface  is  completely  screened  from  the 
incident  pulse.  The  effective  surface  scattering  strength  from  bubbles  near  the 
specular  direction  would  be  estimated  from  Figure  4,  Section  V],  at  about  -30  dB, 
similar  to  backscattering  levels  at  about  the  same  grazing  angles.88  Thus,  the 
scattering  differential  cross  section  per  unit  area  of  the  surface  (see  Section  II)  is 
roughly  a  constant  at  about  ct  -  10-3.  The  fraction  F  of  the  incident  energy 
scattered  into  the  solid  angle  H  about  the  specular  direction  is  given  by 

F  =  a  H/sin  G  .  (B8) 

where  9  is  the  specular  grazing  angle.  If  we  associate  with  the  bounce  path  the 
scattered  energy  within  a  cone  of  half  angle  a  of  the  specular  ray,  then 
0  =  2rr  (1  -  cos  a).  Considering  the  discussion  in  Section  V.A.  it  would  appear 
reasonable  to  choose  the  full  angle  2a  -  2CP,  or  a  =  10°.  However,  the  actual  beam 
widths  and  pulse  length  for  a  particular  application  should  be  taken  into  account  in 
making  such  an  estimate.  The  maximum  bounce  loss  is  then  given  by 
SBLj^^dB)  =  - 10  log  F .  With  the  typical  value  of  6  =  10°,  we  find  SBLfcm,  -  33  dB. 

Alternately,  a  theoretical  estimate  can  be  made  from  the  work  of  McDaniel  and 
Gorman.69  From  Eq.  (15)*  of  Ref.  69  we  have  in  the  limit  of  large  bubble  densities 

o  =  sin  G  ,  (B9) 

8tt  o 


•  A  typographical  error  in  Eq.  (15)  of  Ref.  69  has  been  corrected  (see  Ref  B0),  which  reduces  the  result  by 
7T*  Also,  inclusion  of  Mcdwin's  factor®2  of  1.7  in  attenuation  [Eq.  (10),  Ref.  69]  as  well  as  in  scattering  [Eq 


(8),  Ref 


yields  an  additional  reduction  by  a  factor  of  1.7  Eq.  (39)  can  also  be  written 


1  s 

a  =  -= —  sin  9,  where  S.  and  S  are  the  total  scattering  and  extinction  cross  sections  The  cross  sec- 

07T  O  5  e 


tion  notation  here  follows  Ref  52  and  differs  from  Ref  69 
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where  the  damping  factor  6  depends  weakly  on  the  bubble  radius.08  Note  that  with 
the  use  of  Eq.  (B9)  in  Eq.  (B8),  the  sin  9  factor  cancels.  With  the  use  of  the  typical 
value80  d  =  0.07  and  with  the  same  fi  as  before,  Eq.  (B8)  and  Eq.  (B9)  lead  to 
SBL,^  2  31  dB. 

Thus,  we  adopt  the  estimate  SBl,^  =  30  dB. 
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